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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TECHNICAL MEMORANDUM 1353 


SOME PROBLEMS ON THE THEORY OF CREEP* 
By Y. N. Rabotnov 


Section 1. The term creep of metals is applied to the phenomenon in 
which, at temperatures beyond a certain limit, the metal subjected to a 
load slowly undergoes deformation with time. For the case of steel, the 
creep phenomenon must be taken into account at temperatures above 400° C. 
Very slow deformations for a prolonged period are cumulative and lead 
either to inadmissible changes in the dimensions of a structural part or 
to its failure. 


It is important to note that failure due to creep occurs for very 
small strains considerably less than those in static rupture. 


In the design of steam power units, boilers and turbines, creep is 
a basic factor which determines the choice of the admissible stresses. 
On account of the extreme urgency of the problem, the creep phenomenon 
has claimed the widest attention of metallurgists, physicists, and to a 
lesser extent, technicians. 


At the same time, however, the theory of creep constitutes part of 
the mechanics of dense media and the mechanical formulation of the pro- 
blem may be given as the following: 


A body is subjected to the action of a given system of forces, or 
initial displacements are prescribed on its surface. It is required to 
find the stress distribution in the body and the changes of its defor- 
mations with time. 


Such a statement of the problem immediately raises the following 
question: What tests should be set up in order that the mechanical 
characteristics of creep (certain constants or functions) may be deter- 
mined? Is it sufficient for this purpose to make use of the generally 
accepted methods of testing, or is it necessary to supplement them? 


For the solution of the problem of creep as thus formulated, a 
mechanical theory of creep is required. Such theory, at the present 
state of knowledge of the physics of the process, must necessarily bear 


* "Nekotorye Voprosy Teorii Polzuchesti." Vestnik Moskovskovo 
Universiteta, No. 10, 1948, pp. 81-91. 
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an external, over-all character based, not on the investigation of 
microprocesses, but on the results of mechanical tests. The existing 
physical theories are as yet far from providing a quantitative descrip- 
tion of the process in all its complication. 


The usual method of creep tests is that of obtaining the strains 
for a constant load. For small deformations it may be assumed that to a 
constant load there corresponds a constant stress. 


The results of tests are generally represented in the form of creep 
curves (fig. 1), the time being laid off on the axis of abscissas and the 
strain € on the axis of ordinates. 


The intercept €, represents elastic deformation if the stress o 


does not exceed the elastic limit of the material. Often in place of 
the total strain there is laid off the plastic strain 


p=e- o/E 
Many attempts have been made to give an analytical expression for 


the creep curves. The different equations proposed may be divided into 
two groups: 


l. p = S(o)T(t) | (1.1) 


2. p=e(o)8(t) + S(o)t  (8(*) = 1) (1.2) 
In writing equation (1.2), the essential assumption is that the rate of 
creep tends towards a constant value with time, that is, the creep curve 
has an asymptote. The function g(o) is the intercept ¢' on figure l. 
For the function S(o), the following equations have been proposed: 


S(0) = Aon (Bailey ) 
Sto) = veo! (Ludwik) 
Sto) = y(eS/H s) (Soderberg) 
S(o) = 2v sh = (Nadai) 


Sto). = boe e (Oding) 
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For the function T(t), the best equation is apparently the power equa- 
tion 


T(t) =t™ (0<m<1) 
If equations of type (1.2) are used, the function 0(t) assumes 
the following forms: 


g(t) = 1- e Pt (Mac Vetty) 


o(t) = — (Oding) 


For g(0), it has not as yet been possible to establish a law. 


The series of creep curves for different values of o give the 
representation of the functional relations between the three variables 
o, €, and t. If e and o are laid off on the coordinate axes, a 
series of curves is obtained, which is shown in figure 2 and character- 
ized by different values of t. Figure 2 refers to the test data of 
Robinson which are unique in that these tests extended for 100,000 hours 
(from March 27, 1931 to October 8, 1942). We worked over a very large 
number of the tests in the series by this method and obtained in all 
cases, with an accuracy not exceeding the limits of the experimental 
accuracy, affinely related curves in the oe plane. 


On the basis of this result, the following formula describing the 
law of creep for constant stress is proposed: 
p(e) = [1 + 4(t)] o (1.3) 


For the function G(t), good results are generally given by the following 
expression: 


l-a 


G(t) = + 


where the coefficient a fluctuates about the value 0.7, while the 
coefficient X changes very strongly for different materials. 


The scatter in transferring all the points on to the curve t=0 
(fig. 2) is found to be not greater than the scatter of the modulus of 
elasticity of the specimens on which the given series of tests was con- 
ducted. 
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Section 2. The mechanical theory of creep for the case of a single 6 
axis must establish such relation between € and 0, containing the 
time or time operators, which would permit predicting the course of the 
process varying with time. In particular, a second extreme case of the 
one-dimensional problem is the problem of relaxation. The latter is the 
process of decrease in stress in a rod the length of which remains con- 
stant. Test data on relaxation are very meager, only those data being 
of value which, together with the curves of relaxation, give the curves 


of creep for the same material. The only reliable data in the general o 
literature are those published by Davis in 1943 on copper. W 
The various mechanical theories of creep existing at the present 
time may be divided into three groups: 
(1) Theory of constant rate: Assuming the existence of an asymp- 
totic curve of creep, the curve is replaced by a straight line parallel 
to the asymptote and intercepting the segment €, = E/E on the axis of 
ordinates. Then 
p = S(o) (2.1) $ 


The theory of constant rate assumes this as the true relation for any 
conditions. In particular, for € = constant, p = -5/E, and from equa- 
tion (2.1) there is readily obtained the law of relaxation which grossly 
contradicts test data, since the neglect of the primary creep (the curvi- 
linear part of the curve et) is not permissible in the problem of 
relaxation. In the case of the nonuniform stress state, equation (2.1) 
leads to very great difficulties and even the problem of the pure bending 
of a rod of rectangular cross section is not solvable. 


The greater number of authors employing this theory take still 
another step and neglect the elastic deformation. The fundamental equa- 
tion is then the following: 


é = S(0) (2.2) 


Equation (2.2) is generally put at the basis of the theory of sec- 
ondary creep widely applied in technical computations. Below shall be 
given another basis for this theory which considerably generalizes it. 
For the present we may note that from the point of view of equation (2.2) 
the problem of relaxation has no significance. 


(2) Theory of aging: The theory of aging postulates the existence 
of a definite relation between o, €, and t: > 
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f(o,e, t) = (2.5) 


or between p, 0, t 


F(p, O, t) 


It is easily shown that any theory, the fundamental equation of which 
contains the time explicitly, is contradictory. The physical law must 
be invariant relative to a time origin. Applying the theory of aging 
to the successive loading and unloading may yield absurd results. 


O (2.4) 


The theory of aging in the form of equation (2.3), however, gives 
for smoothly varying loads satisfactory agreement with experiment. The 
methods of computation based on it are relatively simple and at the same 
time permit taking into account all the characteristic experimental 
curves which may be obtained in tests. Hence, one of the variants of 
the aging theory can be recommended as a technical method of computing 
structural parts working under the conditions of high temperatures. This 
point of view has been developed by us in a paper presented at the 
session of the Soviet Academy of Sciences in March 1948. 


(3) Theory of strain hardening: This theory postulates the exist- 
ence of an unvarying relation among the rate of plastic deformation, its 
magnitude, and the magnitude of the stress: 


®(p, p, 0) = 0 (2.5) 


Methods exist for the graphical construction of the relaxation curve 
by a given family of creep curves on the basis of hypothesis (2.5). The 
analytical formulation of this theory, the choice of the functional 
relations for which it is possible to integrate equation (2.5), is found 
to be very difficult. 


Section 3.. The theories of creep enumerated in section 2 are not 
capable of explaining a number of phenomena observed during experiment. 
There is first of all the case for the strain-hardening effect. A 
specimen initially strained by a large force evidences creep to a con- 
siderably less extent than a specimen not initially strained. If the 
specimen was tested under stress 0] and 0; is decreased to do, the 


creep practically vanishes. If, however, the stress do was initially 


imposed, the creep for this stress may be very marked. There are no 
serious investigations of this problem.in the literature. The published 
data undoubtedly give a qualitative account of the phenomenon but the 
quantitative aspect still awaits investigation. 
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The second effect is the so-called reverse creep. If a specimen is è 
subjected to a creep test at constant load and the load is then removed, 
the specimen immediately shortens by the amount of the elastic elongation 
but the process does not stop there. In the course of time the specimen 4 
continues to shorten, returning in this way a part of its residual 
deformation. This phenomenon of the type of elastic aftereffect has been 
subjected to a careful experimental study but not one of the previously 
enumerated theories provides an explanation for it. 


ASIEK 


The theory of creep proposed by us represents an extension of the 
theory of elastic heredity of Volterra to plastic deformation. In the 
same way as the elastic heredity develops about a straight line o = Ee 
in the g€ plane, the plastic heredity, or creep, develops about a 
certain curve in this plane. It will be inconvenient, in what follows, 
to lay off as usual the value of the stress o on the ordinate axis and 
we shall therefore take a certain fictituous plane QE and a curve 
p=«p(€) in this plane (fig. 4). 


The curve (€) represents the ideal curve of strain with the 
exclusion of the time factor (actually never realized). For active 
processes, that is, those accompanied by motion along this curve upward, 
the fundamental law is written in the following manner: 


(e) = (1 + K*) o (3.1) 


where K* is the integral operator of Volterra, that is, 
t 
K*o = K(t - T) o(T) aT 
O 


Here and in what follows, use will be made of the notations and 
results of our previous paper (ref. 1). 


The increase in the stress o is no longer a criterion of the 
activity of the process as in the theory of plasticity but such criterion 
is given by the increase in e. 


For unloading processes, it is necessary in the left-hand side of 


equation (3.1) to introduce in place of (€) the value of the ordinate 
of the linear unloading AB expressed as a function of €: 


E (e - €') +g' = (1 + K*) o (3.2) 
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If T* is the solving operator, 


l » 
z Belial 


there follows from equation (3.1) 


o =(1-T*) oe) 


Setting 


K* + 1=G(t) and T* +1 = R(t) 


the law of creep will then be 


p(e) =[1+4(t)] o, (3.3) 


and the law of relaxation 


o =[1 - R(t)] (e,) (3.4) 


Equation (3.3) entirely agrees with that obtained from tests, for- 
mula (1.3). The expression given in section 1 for G(t) shows that the 
operator K* possesses a singularity of the type of the Abel operator. 
In the simplest case we may assume 


$ 
K* = XP a 


where I*_, is an operator with kernel (t - T)"*/T(1-a). Introducing 


the 3-operators employed in the previously cited paper (ref. 1) yields 
a very close agreement with test results, since the 3-operator contains 
an additional constant. [NACA Reviewer's Note: The 3-operator is an 


e on n(1-+a) 
a B (t -=T 
operator with the kernel 3,(B, t - 1) = (t - 1) ) oe, 
a 5 Filet il +4 


The additional constant referred to is B.] 


If, as is usually the actual case, the body is acted upon by con- 
stant loads or the displacement of its points is maintained by constant 
stresses, it is necessary in actual computation to deal with the two 
functions G(t) and R(t) which may be determined from tests forgetting 
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about their origin from the kernel and as resolvents of the same inte- 
gral equation. 


This variant of the technical theory of creep is a forward step as 
compared with the aging theory which retains equation (3.3) both for 
creep and for relaxation. 


The phenomenon of reverse creep (fig. 5) will now be considered. A 
stress o is applied to the specimen and it undergoes the instantaneous 
deformation o/E (point A); in the course of time T, creep ocgurs along 
the curve AB; at the instant T the stress is removed and at the same 
time there is removed the elastic deformation BC = OA = 0/E. In the 
time elapsed 9 after the unloading there is also removed the defor- 
mation €, (point D). This process will be followed in the plane (pe 


(fig. 6). 


The instantaneously applied load corresponds to the motion along 
the curve from point O to point A; during time T the strain increases 
and the function ple) up to point B; the creep process is described 
by equation (3.1). At point B 


(eg) =L1 + G(T)] o (3.5) 


With instantaneous unloading we drop to point C on the segment 
BC = AO, the reverse creep corresponding to the motion along this line 
up to the point Dat instant 8. [NACA Reviewer's Note: It is also of 
interest to examine the effect of an instantaneous increase in load 
occurring after a period of constant-stress creep. If €, is the strain 


at the end of the creep period and 6€ is an instantaneously imposed 
strain increment, then, according to equation (3.1), there is an 
instantaneous stress increment given by 


50 = p(co + 5e) - OE) (a) 


This result may be at variance with the facts, for recent experimental 
evidence, including tests on the propagation of plastic waves in bars 


subject to creep, suggests that after some creep has taken place materials 


behave elastically for small instantaneous increments of stress, that is, 


Do = E SE (b) 
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Only for sufficiently small €g would equation (a) reduce to equa- 
tion (b).] Here equation (3.2) must be used: 
E (€ - cp) + Pp = [G(T + 0) - G(@)] o (3.6) 


[NACA Reviewer's Note: The coefficient o in equation (3.6) was 
incorrectly omitted in the original.] 


It is seen from figure 6 that 


Making use of equations (3.5) and (3.6) yields 


4 = g {o(2) + a(o) - G(T + 0) (3.7) 


The linear dependence, of the reverse creep on the stress and the 
symmetry of its dependence ` on T and 6 is well confirmed by experi- 
ment (ref. 2). 


op portant particular case of the relation (3.1) is obtained when 
ple) = e, and K* > q* The law of creep is then 


BE X = a 
Sa ie = Tien a) 9 O 


For large t, the second term in the brackets is the dominating one and 


e = Aot” (3.8) 
where 
B? B 


If 1-a=8, formula (3.8) gives, in the limit, a constant creep 
rate. 


Section 4. The application of the theory presented in the preceding 
section to the problem of pure bending is considered herein. For sim- 
plicity, we restrict ourselves to the case of a rod of rectangular cross 


section, since we are concerned with the theoretical side of the problem. 


Assuming the hypothesis of plane sections, set 
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e = z/p 


where p is the radius of curvature of the bent axis, and z the 
coordinate measured from the neutral axis and varying from -h to +h. 


By equation (3.1) 


(1 +K*) o = de (4.1) 


Multiplying by bzdz, where b is the width of the section, we 
integrate from z =0 to z = h, and multiply the result by 2 (the 
function ple ) is analytically determined only for positive values of 
the argument; it is continued as an odd function in the region of nega- 
tive values). Since the time and space operators are interchangeable, 


h 
(1 + K*) M = 2b (2) zăz 


O. 


where M is the bending moment. Introducing the notation 


x 
ii 
x2 p(x) xax = p(x) 
i O 
yields 
(1 + K*) M = 2bh* 19 (4.2) 


The graph for the function o(h/p) can be easily constructed. It 
is thus always possible to find the magnitude of the curvature from the 
given moment M for a given instant of time and then, by solving inte- 
gral equation (4.1), to determine the stress distribution. 


If the function (e) is a power function, the problem is consi- 
derably simplified. Let P(x) = axP. Then 


P(x) = a xP 


Dividing equation (4.1) by (4.2) yields 


S92 
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whence 


Ree (zy M (4.3) 


The distribution of the stresses is found to be independent of the time. 
We thus have in a certain sense a secondary creep although the rate of 
deformation is not constant as is clear from equation (4.2). 


The result, equation (4.3), was obtained by Davis (ref. 3) from 
entirely different considerations. He made use of the theory of strain 
hardening, using equation (2.4) in the form 


p = Bp'o 


As he found difficulty in accurately solving the problem as thus 
posed, Davis identified the plastic deformation with the total defor- 
mation and arrived at a creep law of the type (3.8) from which he readily 
obtained equation (4.3) in the case of bending. The experimental verifi- 
cation conducted by Davis well confirms the theoretical result. 


Section 5. The generalization of any of the theories of creep to 
the three-dimensional case may be effected if there are taken as a basis 
the equations of the theory of small elasto-plastic deformations: 


40. 
i 


Ox = O = 73 (5.5) 


The change in volume at high temperatures may be neglected with 
greater justification than for normal temperatures since, starting from 
4000, the Poisson coefficient for steel assumes a value of the order of 
0.45 to 0.47. 


The intensity of the stresses and the strains which are determined, 
following A. A. Ilyushin, as 
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2 2 
0; =% (0, - ene + (oy - 07) + (0, - 07) +6 (Txy + Tyz + Taa) 


£G9Z 


are assumed to be connected with each other by the same relations as the 
stress and strain in the one-dimensional problem. 


By the theory herein, 


plei) = (1+ K*) 0 (5.2) 


The problem of a pipe under the action of an internal pressure is 
very simply solved on the assumption of the absence of axial deformation. 
Let 0, and Og be the radial and transverse stresses, tr = du/dr 


and €g = u/r the corresponding strains expressed in terms of the radial 


displacement u, b and a the external and the internal radius of the 
pipe, and q the internal applied pressure. Since €, = 0, there 


follows from the condition of incompressibility 


du u 
dr te 0 
whence 
e _ €n =V¥2 £8 € _ eat 
BESS E ee a 


where e is an as yet undetermined funetion of time. Equations (5.1) 
give [NACA Reviewer's Note: In the following three equations, an error 
in sign appearing in the original has been corrected. |] 
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whence 


2 
O - O = = (O. - Deo 
r O NE L ( ) 


Substituting (5.3) in the equation of equilibrium 


do, Oi F 03 
de "e 0 
yields 
do o 
r 2 1 | 
EI WEJ (5.4). 
Multiplying by (1+ K”); integrating from r =a, and noting that 
+ ea? 
(1 + K*) o; =(e,) (E 
T 
yield 
T 
. ae ea”) ar 
SES o (5): 
Setting 
x 
dx 
p(x) = = (x) 
O 
gives 
(1 + K*)(o. + a) = lo(e) - (e a (5.5) 
NE r 
For r = b, 0, = 0. From this condition, the integral equation for 
determining the function of time e is obtained: 
| 1 al | 
(1+K*)q = = |$(e) - Qe 5 (5.6) 
NE p? 
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The case where (x) = axP here too leads to considerable simplification; 
in this case 


(x) =5 xP 


[NACA Reviewer's Note: The parameter a in ax® is evidently being 
taken as l. 


cgiz 


From equation (5.5), 


ap 
(1 + K%)(0, + q) EF ep -S 


(1 + K*)q = 


+ B a2P 

e l -= E 

B /3 p2P 
Dividing one equation by the other cancels the time function e on 


the right side. The integral operator (l + K*) likewise cancels and a 
time-independent law for the stress distribution is obtained: 


2B _2B 2B 
a r - b 
> ae) „28 _ ¿28 (5.7) 


The velocities are not constant and they are found from equation (5.6) 
where, if q is given, the case reduces to quadratures and the succes- 
sive determination of the function e from the graph. If e is given, 
the problem reduces to the solution of an integral equation for which it 
is necessary to know the resolvent of the kernel. 


Section 6. Proceeding to the general case of the three-dimensional 
problem, we restrict ourselves to the consideration of those processes 
in which the stresses vary simultaneously in proportion to a parameter 
x(t), while the strains vary in proportion to a parameter u(t). Then 


0% = No, 4.4 e E, = HE x 
where 6%; o.. ES . . . are functions only of the coordinates. 
[NACA Reviewer's Note: The bars which appear above 0, and €, in this 


clause were incorrectly omitted in the original.] Equations (5.1) give 


Moy - 0) = 3 03 = 
i 
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The stresses 0, are connected with the strains by the relation (552: 
If ©Q(€;) is a power function 
p(€;) = aeh = auf ER (6.2) 


Multiplying both sides of equation (6.1) by 1 + K* gives 


z 
ANA (6.3) 
cd 


[NACA Reviewer's Note: Errors in the subscripts on the right side of 
equation (6.3) appearing in the original have been corrected. | 


Equations (6.3), of which only the first is written out, correspond to 
the integral equation 
(1 + K*) à= pf (6.4) 


and to a system of relations entirely agreeing with the equations of the 
theory of small elasto-plastic deformation: 


EE OE 1 Si — | 
Oy =- 0 = Z a x? Txy z = E xy (6.5) 
1 1 
where 
s, % . zał = | | 
0j =a ej =pl€y) (6.6) 


Systens (6.5) and (6.6) are systems of equations of the theory of 
secondary creep taken in a more general sense than the ordinary. The 
rates are now no longer constant but variable since the factors A and 


u which depend on the time are obtained from the integral equation (6.4). 
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